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Exploring the recent expansion history of the universe promises insights into the cosmological 
model, the nature of dark energy, and potentially clues to high energy physics theories and gravi- 
tation. We examine the extent to which precision distance-redshift observations can map out the 
history, including the acceleration-deceleration transition, and the components and equations of 
state of the energy density. We consider the ability to distinguish between various dynamical scalar 
field models for the dark energy, as well as higher dimension and alternate gravity theories. Finally, 
we present a new, advantageous parametrization for the study of dark energy. 
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I. INTRODUCTION 

The quest to explore the expansion history of the uni- 
verse has carried cosmology well beyond "determining 
two numbers" - the present dimensionless density of mat- 
ter fi m and the present deceleration parameter g 0- 
Observations have advanced so that now cosmologists 
seek to reconstruct the entire function a(t) representing 
the expansion history of the universe. 

A myriad of cosmological observational tests can probe 
the function a(t) more fully, over much of the age of the 
universe (see Sandage 0, Linder [|[ Q], Tegmark ||). 
This paper concentrates on the most advanced method, 
the magnitude-redshift relation of Type la supernovae. 

Just as understanding the thermal history of the early 
universe has taught us an enormous amount about both 
cosmology and particle physics (see, e.g., ||), the recent 
expansion history of the universe promises similarly fer- 
tile ground with the discovery of the current acceleration 
of the expansion of the universe. This involves concepts 
of the late time role of high energy field theories in the 
form of possible quintessence, scalar-tensor gravitation, 
higher dimension theories, brane worlds, etc. as well as 
the fate of the universe. 



II. MAPPING THE EXPANSION HISTORY 

Type la supernovae, or any standardizable candles 
(sources with known luminosity), are excellently suited 
to map the expansion history a(t) since there exists a 
direct relation between the observed distance-redshift re- 
lation d(z) = (1 + z)r)(z) and the theoretical a(t). For a 
flat universe, assumed throughout, 
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where rj is the conformal time, dr) = dt/a, the Hubble 
parameter is H = a/ a and a e = 1/(1 + z) is the scale 
factor at the time of emission, i.e. when the supernova 
exploded. 

The Hubble parameter in general relativity comes from 
the Friedmann equation 



and the conservation condition of each component is 

p/p=-3H(l+p/p) = -3H[l + w(z)}, (3) 

where the energy density is p, the pressure p, and the 
equation of state (EOS) of each component is defined by 
w = p/p. Ordinary nonrelativistic matter has w — 0; a 
cosmological constant has w = — 1. We explicitly allow 
the possibility that w evolves. 

The expansion history a(t) is given by 

t(o)= / da'/(a'H)= dz'/[(l + z')H(z')}. (4) 



The fly in the ointment is that the measured distance 
d(z) is related to, but is not, the desired history rela- 
tion a(t). To translate d(z) into a(t) directly involves a 
derivative of d(z), so noisy data can introduce difficul- 
ties [0, Instead one finds H(z) through the cosmology 
parameters p and w(z). 

Observational evidence for accelerated expansion in- 
forms us that there must be a component with a strongly 
negative EOS - "dark energy" - in addition to matter. 
Then combining equations (Q)-(||): 
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where Cl m is the dimensionless matter density 
87rp m /(3ifQ) and Hq = H(0) is the Hubble constant. 
Eq. (||) can then be used to obtain a(t). 

To find w(z) one could solve the scalar field equation 
for a particular theory but this does not allow a model in- 
dependent parameter space in which to compare models. 
Instead, for generality of treatment, various parametriza- 
tions of w(z) are used. 



A. Linear w(z) 

The conventional first order expansion to the EOS, in- 
cluding the critical property of time variation in the EOS, 
is w(z) = wq + w\z. In this case the exponential in 



(2) resolves to (1 + z) 
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However this grows 
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increasingly unsuitable at redshifts z > 1. Analyzing 
CMB constraints on the distance to the last scattering 
surface at zi ss — 1100 would be problematic. 



B. A new parametrization of the equation of state 

To extend parametrization of dark energy to redshifts 
z > 1, I suggest a new model: 



w(a) = w + W a (l — a) 

= w + w a z/(l + z). 



Here the 
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resolves to 



This new parametrization has several advantages: 1) 
a manageable 2-dimensional phase space, 2) reduction to 
the old linear redshift behavior at low redshift, 3) well 
behaved, bounded behavior for high redshift, 4) high 
accuracy in reconstructing many scalar field equations 
of state and the resulting distance-redshift relations, 5) 
good sensitivity to observational data, 6) simple physical 
interpretation. 

Beyond the bounded behavior, the new parametriza- 
tion is also more accurate than the old one. For example, 
in comparison to the exact solution for the supergrav- 
ity inspired SUGRA model || it is accurate in matching 
w(z) to -2%, 3% at z = 0.5, 1.7 vs. 6%, -27% for the linear 
z approximation (the constants w\, w a are here chosen 
to fit at z = 1). Most remarkably, it reconstructs the 
distance-redshift behavior of the SUGRA model to 0.2% 
over the entire range out to the last scattering surface 
w 1100). 

Note that dw/d\n(l + z)|z=i = w a /2; one might con- 
sider this quantity a natural measure of time variation 
(it is directly related to the scalar field potential slow 
roll factor V' /V) and z = 1 a region where the scalar 
field is most likely to be evolving as the epoch of matter 
domination changes over to dark energy domination. 

The future Supernova/ Acceleration Probe (SNAP: 
pof) will be able to determine w a to better than ±0.55 
(one expects roughly w a ~ 2u>i), with use of a prior on 
n m of 0.03, or to better than 0.3 on incorporating data 
from the Planck CMB experiment |ll| . For the advan- 
tages of combining supernova and CMB data see Jl2| . 
The CMB information can be folded in naturally in this 
parametrization, without imposing artificial cutoffs or lo- 
cally approximating the likelihood surface as required in 
the w\ parametrization. In fact, the new parametrization 
is even more promising since the sensitivity of the SNAP 
determinations increases for wq more positive than —1 or 
for positive w a (see, e.g., the values quoted above 
were for a fiducial cosmological constant model. For ex- 
ample, SUGRA predicts w a = 0.58 and SNAP would put 
error bars of a(w a ) ps 0.25 on that; this would demon- 
strate time variation of the EOS at the 95% confidence 
level. Incorporation of a Planck prior can improve this to 
the «99% confidence level: a(dw/d\n(l + z)| z =i) « 0.1. 



C. Expansion and density histories 

Figure |l| shows the mapping of the expansion history 
of various models and the constraints that SNAP data 
would impose. The models include both dark energy and 
alternative explanations for the acceleration discussed in 
Sec. HI. Regarding the data constraints, it is important 
to note the presence of correlation between the param- 
eters f2 m , wo, and w\ or w a (see, e.g., [fn|). Despite 
SNAP being able to determine each individually to high 
precision, e.g. f2 m to 0.03, w to 0.05, w\ to 0.3 (each 
marginalized over others), the correlations among them 
relax the tightness of the constraint SNAP would place 
on t he ex pansion history. This is unavoidable (but see 
Sec.|ilD|). 
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FIG. 1: Mapping the expansion history through the super- 
nova magnitude-redshift relation can distinguish the dark en- 
ergy explanation for the accelerating universe from alternate 
theories of gravitation, high energy physics, or higher dimen- 
sions. All three models take an Am = 0.3, flat universe but 
differ on the form of the Friedmann expansion equation. 



Figure g gives the equivalent density history. Note 
that the matter dominated epoch is characterized by 
p(z) ~ (1 + z) 3 and so has a slope of 3 in this plane. 
The deviation from this due to the recent epoch of dark 
energy domination can clearly be seen (cf. Tegmark 0). 



D. Conformal time history 

One method of incorporating the advantages of both 
mapping approaches - the generality of parametrization 
and the directness of reconstruction - is to consider the 
conformal time a{rj). From d = (1 + z)r\ one sees that 
one requires no foreknowledge or local approximation to 
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obtain the scale factor-conformal time relation. The error 
estimation from the observed magnitudes m is simply 
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(8) 



A 1% distance measurement error (cr m — 0.02) given by 
SNAP's limiting systematics becomes a 1% error in a{rj). 

Another virtue of the conformal time history is its 
physical interpretation. Consider the logarithmic deriva- 
tive of the 77(a) curve. This is 



dr) dt _ x 

— — = a— = (a . 
a In a a da 



(9) 



Two simple physical interpretations apply: 1) this gives 
the proper time evolution of the scale factor, 2) this rep- 
resents the conformal horizon scale (aH)^ 1 . Both are 
familiar from analysis of inflationary cosmology - they 
determine when the expansion of the universe enters an 
accelerating phase. 

Figure [| maps the conformal horizon. When this curve 
has positive slope, the expansion is decelerating, a < 
and comoving wavelengths (e.g. of density perturbations) 
enter the horizon. More recently the curve has negative 
slope; this is the signature of accelerating expansion. In 
terms of inflation we would say that wavelengths exit the 
horizon then (the horizon appears to shrink). The dashed 
lines show the reconstruction possible by SNAP observ- 
ing supernovae out to z = 1.7. SNAP clearly has the 
ability to probe not only the current accelerating epoch 
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FIG. 3: The conformal horizon scale characterizes the dy- 
namics. The negative slope part of the curve allows comoving 
wavelengths to expand outside the horizon, or alternately rep- 
resents aH = a increasing, i.e. a > - the signature of infla- 
tion or acceleration. The dashed blue lines show that SNAP 
will map the accelerating phase, the transition, and into the 
matter dominated, decelerating phase of the past universe. 



("late time inflation") but to reach back into the decel- 
erating phase and to map the transition between them. 
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FIG. 2: Mapping the density history with SNAP can clearly 
show both the current accelerating phase and the transition 
to the matter dominated, decelerating epoch. At high red- 
shifts all models obey the matter dominated slope 3 law, but 
deviation from this - the sign of dark energy or alternative 
gravity - is clearly visible for z < 2. 



III. BEYOND DARK ENERGY 

Mapping the physical time evolution of the scale fac- 
tor relies on translating the observations into the Hubble 
parameter H(z). The Friedmann equation (ph related 
this to the matter and dark energy densities. But ide- 
ally we would like to use the data to test the Friedmann 
equations of general relativity or alternate explanations 
for the acceleration besides dark energy. Here we lay the 
framework for such options and give some examples. 



A. Higher dimension theories 

In the braneworld scenario of Deffayet et al. |l4j, grav- 
itation leaks from the 4-dimensional brane we experience 
out into the 5-dimcnsional bulk. The new expansion evo- 
lution equation is 
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where r c = Mj,,/(2Mf) is the crossover scale on which 
gravity leaks into the bulk, defined in terms of the 
usual Planck scale Mp; and the 5-dimensional Planck 
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scale M5. It is convenient to introduce a dimension- 
less energy density $l. rc — 1/ (AH^r 2 ). Flatness imposes 
n rc = (1 - ft m ) 2 /4. 

SNAP data will be able to constrain the parameter 
uncertainties to cr(S! m ) = 0.008, or a(il rc ) = 0.003 for a 
fiducial Q m = 0.3. This corresponds to a determination 
of the crossover scale of r c — 1.43 ± 0.015 Hq 1 . If other 
observations, analyzed within this braneworld picture, 
disagreed with this narrow range then this theory could 
be ruled out. Fig. |l| illustrates how the expansion history 
for a braneworld model differs from that for dark energy. 

Note that other higher dimension models discussed, 
such as the Randall-Sundrum type 2 model [[t5| [l6| and 
so-called Cardassian expansion jl7|] , reduce to the already 
discussed dark energy prescription for the recent expan- 
sion history. 

B. Chaplygin gas 

A very different possibility to explain the accelera- 
tion of the universe is the Chaplygin gas fis] , proposed 
to unify dark energy and dark matter. Its pressure 
p ~ — 1/p gives a solution that at early times behaves 
like nonrelativistic matter and at late times like a cos- 
mological constant. The expansion equation becomes 

[H(z)/H ] 2 = n m (l+z) 3 +{l-n m )^A + (1 - A)(l + zf. 

The factor A is interpreted as the sound velocity squared 
of the Chaplygin gas. As A — > 1, it reduces to the cos- 
mological constant. 

SNAP data will be able to constrain the parameter 
uncertainties to u(Om)= +0.005,-0.017 and a (A) — 
—0.005 for a fiducial fl m — 0.3, A = 1 (though the re- 
sults are fairly insensitive to A) . Note dark matter in the 
form of "quintessence clumps" would require A signifi- 



cantly smaller than unity. 



IV. CONCLUSION 

The geometry, dynamics, and composition of the uni- 
verse are intertwined through the theory of gravitation 
governing the expansion of the universe. By precision 
mapping of the recent expansion history we can hope 
to learn about all of these. The brightest hope for this 
in the near future is the next generation of distance- 
redshift measurements through Type la supernovae that 
will reach out to z w 1.7. 

Just as the thermal history of the early universe taught 
us much about cosmology, astrophysics, and particle 
physics, so does the recent expansion history have the 
potential to greatly extend our physical understanding. 
With the new parametrization of dark energy suggested 
here, one can study the effects of a time varying equa- 
tion of state component back to the decoupling epoch 
of the cosmic microwave background radiation. But even 
beyond dark energy, exploring the expansion history pro- 
vides us cosmological information in a model independent 
way, allowing us to examine many new physical ideas. 
From two numbers we have progressed to mapping the 
entire dynamical function a(t), to the brink of a deeper 
understanding of the dynamics and fate of the universe. 



Acknowledgments 

This work was supported at LBL by the Director, 
Office of Science, DOE under DE-AC03-76SF00098. I 
thank Alex Kim, Saul Perlmutter, and George Smoot 
for stimulating discussions and Ramon Miquel and Nick 
Mostek for stimulating computations. 



[1] 
[2] 
[3] 

[5] 
[6] 

[7] 

[8] 

[9] 



A. Sandage, Ap.J. 133, 355 (1961) 

A. Sandage, Ann. Rev. A&A 26, 561 (1988) 

E.V. Linder, A&A 206, 175 (1988) 

E.V. Linder, First Principles of Cosmology (Addison- 
Wesley, London, 1997) 



M. Tegmark, Science 296, 1427 (2002), |astro-ph/0207199 
E.W. Kolb and M.S. Turner, The Early Universe 
(Addison- Wesley. New York . 1990) 
P. Astier, istro- ph/0008306; 

D. Hut erer and M.S. Tur ner. Phys. Rev. D 64, 123527 
(2001), |astro-ph/001251C 

P. Brax and J. M artin, Phys. Lett. B 468, 40 (1999), 



astro-ph/990504C 



[10] SNAP (http://snap.lbl.gov 



parameter error estimates 
are for 2000 SNe between z = 0.1 - 1.7 and 300 
SNe at z < 0.1 from the Nearby Supernova Fac- 



R. Miquel and N. Mostek; see Miquel, Mostek, & Lin- 
der 2002, in preparation, for a description of the Monte 
Carlo fitter. 



tory (|http: / / snfactory.lbl.goy ) , plus a prior on SI m of 



[11] 
[12] 

[13] 

[14] 

[15] 

[16] 

[17] 
J18I 



Planck: http:/ /astro. estec.esa.nl/Planck 

J. A. Frieman, D. Huterer, E. V. Linder, and M .S. Turner, 

submitted to Phys. Rev. D, istro- ph/020810(] 

J. Well er and A. Albrec ht, Phys. Rev. D 65, 103512 

(2002), |astro-ph/010607S 

C. Deffayet, G. Dvali, G. Gabad adze, Phys. Rev. D 65, 
044023 (2002), |astro-ph/010506^ 

L. Randa ll and R. Sun drum, Phys. Lett. B 83, 4690 



(1999), |hep-th/9906064 

P. Binetruy, C. Deffayet, U. Ellwanger, D. Langlois 

Phys. Lett. B 477, 285 ( 2000). |hep-th/9910219| 

K. Freese and M. Lewis, |astro-ph7o"201229| 

J.C. Fabris, S.V.B. Gongalves, and P.E. de Souza, astro- 



ph/0207430 



0.03. Some numerical results were kindly provided by 



